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Technical Field 

The present invention is related to the field of communication in general, 
and wireless communication, in particular. 

BACKGROUND 

5 Advances in integrated circuit, microprocessor, networking, 

telecommunication and other related technologies have led to wide spread 
adoption of wireless communication, e.g. mobile wireless "cell" phones. In the 
case of wireless communication, such as mobile wireless "cell" phones, typically, 
a mobile wireless "cell" phone (also referred to as a mobile handset), 

10 communicates with a "nearby" service station (also referred to as a base station), 
which relays the communication signals for the mobile handset. The 
service/base station provides the relay service for all mobile handset in its 
coverage area (its "cell"). Thus, a service/base station typically receives and 
processes communication signals from a number of mobile handsets in its "cell" 

15 concurrently. 

Use of multiple antennas at the service/base station for receiving the 
communication signals from the mobile handsets have become popular, as it has 
several advantages in terms of enhancing the capacity and throughput of the 
wireless communication system. Various signal processing techniques are 

20 employed to process the received communication signals, including but not 
limited to "space-time" processing techniques. 

Among the space-time processing techniques, beamforming is one of the 
promising areas of interest for enhancing the strength of signals received from a 
desired direction. One known technique is the employment of a known training 

25 sequence to estimate the optimum weights (e.g. using least mean square 

(LMS)), for beamforming to a desirable direction. Other known techniques for 
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estimating the directions of arrival (DOAs) include employment of the Bartlett 
processor or the MUSIC (Multiple S/gnal Classification) technique. 

Training has the disadvantage of incurring overhead in the throughput of 
the system, and convergence may take longer time than the time available to 
5 make the determination. The latter techniques require a large number of 
snapshots of the received signals to provide a good estimate of a correlation 
matrix reflective of the correlation (or the lack thereof) of the received signals 
from the independent signal sources. 

Additionally, while each communication signal typically has a number of 
10 multipaths, due to environmental factors, such as reflection off structures and so 
forth, these techniques typically estimate the DOA based only on the most 
dominant multipath of a signal. 
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BRIEF DESCRIPTION OF THE DRAWINGS 
Embodiments of the present invention will be described by way of the 
accompanying drawings in which like references denote similar elements, and in 
5 which: 

Figure 1 illustrates a communication environment suitable for practicing 
the present invention, in accordance with one embodiment; 

Figure 2 illustrates a portion of the operational flow for determining the 
directions of arrival of J signals wirelessly transmitted from J signal sources, in 
1 0 accordance with one embodiment; 

Figure 3 illustrates a portion of the operational flow for determining the 
directions of arrival of L multipaths of J signals wirelessly transmitted from J 
signal sources, in accordance with one embodiment; 

Figure 4 illustrates a computer system suitable for use to practice one or 
1 5 more aspects of one or more of the signal processing methods of Fig. 2-3, in 
accordance with one embodiment. 
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DETAILED DESCRIPTION OF EMBODIMENTS OF THE INVENTION 

Embodiments of the present invention include but are not limited to 
5 methods and apparatuses for determining the direction of arrivals of a number of 
signals wirelessly transmitted. Embodiments of the present invention also 
include methods and apparatuses for determining the direction of arrivals of 
multipaths of a number of signals wirelessly transmitted. 

In the following description, various aspects of embodiments of the 
10 present invention will be described. However, it will be apparent to those skilled 
in the art that other embodiments may be practiced with only some or all of the 
described aspects. For purposes of explanation, specific numbers, materials and 
configurations are set forth in order to provide a thorough understanding of the 
embodiments. However, it will be apparent to one skilled in the art that other 
1 5 embodiments may be practiced without the specific details. In other instances, 
well-known features are omitted or simplified in order not to obscure the 
description. 

Various operations will be described as multiple discrete operations in 
turn, in a manner that is most helpful in understanding the embodiments, 
20 however, the order of description should not be construed as to imply that these 
operations are necessarily order dependent. In particular, these operations need 
not be performed in the order of presentation. 

The phrase "in one embodiment" is used repeatedly. The phrase 
generally does not refer to the same embodiment, however, it may. The terms 
25 "comprising", "having" and "including" are synonymous, unless the context 
dictates otherwise. 
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Referring now to Fig. 1 wherein an overview of a communication 
environment suitable for practicing the present invention, in accordance with one 
embodiment, is shown. As illustrated, for the embodiment, communication 
5 environment 100 includes J mobile handsets 102a-102j and base station 106, 
communicatively coupled to each other. Base station 106 relays communication 
signals for mobile handsets 102a-102j including receiving signals wirelessly 
transmitted by the handsets. The received signals, among other things, may be 
processed and forwarded to another signal processing node (not shown) of a 

10 wireless communication service network, to which base station 106 is coupled, 
and/or to a signal processing node (not shown) of a PSTN (Public Switched 
Telephone Network). 

For the embodiment, base station 106 includes N antennas 108a-108n, 
RF unit 110, and signal processing unit 110, coupled to each other as shown (RF 

15 = Radio Frequency). Antennas 108a-108n are employed to transmit to, and 
receive signals from mobile handsets 102a-102j. Additionally, antennas 108a- 
108n may be employed for other purposes. Antennas 108a-108n may also be 
referred to as sensors. For the purpose of this application, the two terms are 
synonymous. 

20 RF unit 110 is employed to down convert RF signals received by antennas 

108a-108n into baseband signals, or up convert baseband signals into RF 
signals for transmission by antennas 108a-108n. 

Signal processing unit 110 is employed to process the down converted 
baseband signals, and process outbound signals for up conversion. For the 

25 embodiment, signal processing unit 110, includes in particular DOA Estimation 
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unit 112 and Beamforming unit 114, coupled to each other and to RF unit 110 as 
shown. 

In various embodiments, DOA Estimation unit 112 estimates the DOAs of 
the J signals (J being an integer), to be described more fully below referencing 
5 Fig. 2. For these embodiments, Beamforming unit 1 14 forms the corresponding 
weighted output signals based at least in part on the DOAs of the J signals 
estimated by DOA Estimation unit 112. 

In other embodiments, DOA Estimation unit 112 estimates the DOAs of 
the L multipaths (L being an integer) of the J signals, to be described more fully 
10 referencing Fig. 3. Some embodiments support one of the two methods to be 
described referencing Fig, 2 and 3. Other embodiments support both. For those 
embodiments supporting both methods illustrated by Fig. 2-3, the embodiments 
may also support configuration of one of the two methods as the current 
operating method. The configuration may be supported in a static or dynamic 
15 approach. 

For the latter embodiments, Beamforming unit 114 forms the 
corresponding output signals based at least in part on the DOAs of the L 
multipaths of the J signals estimated by DOA Estimation unit 112. More 
specifically, Beamforming unit 114 determines a number of corresponding 
20 weights and forms the corresponding output signals based at least in part on 
weighted contributions of the L multipaths of the J signals estimated by DOA 
Estimation unit 112. 

Except for the advantageous manners base station 106 acquires signals 
104a-104j, mobile handsets 102a-102j and RF units 110 represent a broad 
25 range of these elements. A computer system suitable for hosting a software 
implementation of DOA Estimation unit 112 and Beamforming unit 114 will be 
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further described below referencing Fig. 4. The software implementation may be 
developed employing any one or more of a number of programming languages. 
However, the present invention anticipates all or portions of DOA Estimation unit 
112 and Beamforming unit 114 may be implemented in hardware, using e.g. one 
5 or more Application Specific Integrated Circuit (ASIC) or Field Programmable 
Logic Devices. One of ordinary skill in the art would be able to do so, based on 
the description provided herein. 

For the remaining descriptions to follow, and for the claims, the following 
10 conventions are employed: 

- boldface-capital letters, such as A, M, represent matrices or 

subspaces, 

- boldface-small letters, such as v, represent vectors, and 

- non-boldface letters, such as "s", represent scalars. 

15 

Referring now to Fig. 2, a flow diagram illustrating the operation flow of 

DOA Estimation unit 112 for determining the DOAs of the J signals in accordance 

with one embodiment, is shown. The embodiment assumes the followings: 

The J independent signals impinge on the N antenna array in J distinct 
20 directions 0, f ...0 y , where the angles 6 } are measured with respect to the endfire 

direction. 

The output signal vector for a single snapshot of the received signal is 

given by 

(1) x(f) = As(f) + n(0 

25 where x(t) is the N x 1 column vector written as, 

(2) x(o=Mo,..,*i(or 
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(3) A = [a(0, ),...a(0 7 )], where A is the NxJ directional response 

matrix 

(4) a(0,)=|l . . e- iiL ' lWam0 'f 9 j = l 9 ... 9 J 

are the signal directional vectors. Each a(0 y ) vector has the 

5 dimension Nxl . 

(5) s(0 = [s t (t) 9 ... 9 Sj(t)Y , s(0 is the 7x1 vector that contains the 

/signals transmitted by J independent sources. 

k = 27f 0 /c 9 f 0 being the center frequency, and c being speed of 

electromagnetic wave, k being the wavelength, and 
10 d is the spacing between antenna elements and t corresponds to the 

time index.. 

Thus, x(0 is the linear combination of J-linearly independent 
Vandermonde vectors, each corresponding to a source direction, and x(t) may 
be considered as 

1 5 (6) x(0 = a(0, )s t (0 + . . . + a(0 y )s y (0 

In Equation (3), matrix A represents a basis that spans a J-dimensional 
subspace, to which x(/) belongs. Accordingly, x(/)in its decomposition has one 

of the directional vectors as its element. 

Hence, there exists a correlation such that the intersection of the one- 
20 dimensional subspace spanned by x(t) (output signal vector for a single 

snapshot) and any of it's J-signal directional vectors is non-zero. 

Therefore, if 0,. does not correspond to one of true source directions 

i.e.0. g(0 iv ..,0 y ), the intersection between x(0and a(0 / )will be zero. 

Accordingly, the computational problem is formulated as 

25 (7) x(f) n a(0. ) =0, for 0. £ (0, v ..,0 y ) 

(8) x(0 n a(0 f ) ?t 0 for 0. e (6, 0 y ) 
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The Nx2-matrix formulated as D(/) = [x(*),a(0,)] will have rank deficiency 
or close to rank deficiency if 0. € (0, ,.„,0 y ) . 

Hence, in determining the DOA of the J signals 0,,...0 y , the process of Fig. 
2 first selects a trial direction for a signal source y, - 0. , block 202. 

On selecting a trial direction, the process proceeds to compute 



On computing the first coefficient and the first orthonormal vector, the 
process proceeds to compute 



determines whether the evaluation yields a peak value for the function. If 
so, the trial direction is considered to be one of the DOAs of the J signals. 

At block 210, the process determines whether additional directions are to 
be determined. If so, the process returns to block 202, and continues from there 
as earlier described. If not, the process terminates. 

Referring now to Fig. 3, a flow diagram illustrating the operation flow of 
DOA Estimation unit 112 for determining the DOA of L coherent multipaths of J 
signals in accordance with one embodiment, is shown. The embodiment 
includes an initial determination of the DOAs of the J signals, and the determined 



a first coefficient (9a) r u = \\x(t)\\ , and 

a first orthonormal vector (9b) q, , block 204. 



a second coefficient (10a) r n = qf a(~ 0.) , 

a second orthonormal vector (10b) q 2 = a(~ 0,) -r 12 q, , and 
a third coefficient (10c) r 22 = q"a(~ 0.)| , block 




Then, the process proceeds to evaluating a function 



(11) 5(0) = — for ~0., and 

r 22 
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DOA of the J signals are then employed, along with other information, to 
determine the DOAs of the L multipaths of the J signals. 

The signal corresponding to /th source and L coherent multipaths at the 

base station is given by 

5 (12) Xj </) = H ™ e »- fTl,) s{t - 1, ) 

Extending the above equation for N-element antenna array, the following 
relationship is obtained. 

(13) x . (o = ^ vtf, yR,*"****,-*, > s(t _ Tji ) 

5 

-\27iCose }l HnCosOji i27r(N-\)Cos0ji 

y (0 jl ) = [l e " x e x e 1 ] T is the 

where 

10 array response vector for jl-th multipath component. The notations in the above 

have the following meanings. 

R yl in the signal strength or amplitude of the /'th multipath. 

fjCosOj! is the Doppler shift of the ;7th multipath , 

Tjl = Vjl /c is the time delay of the of the jl-th multipath, 

15 r n is the range of the yY-th multipath, and 

c is the speed of electromagnetic wave, 
d fl is DOA corresponding to jl-th multipath. 

L is the number of dominant multipaths for j -th source. 

Under the assumption that multipath delay spread 

20 T — max(r jf ) - min(x Jt ) « 1/B , 

where B is the bandwidth of the signal, 
Sj(t) , the base band signal of the jth source can be modeled as a 

narrowband signal, that is, s(t - x j{ ) » s(t - r 0 ) , where t 0 g [min(x ^maxix ^)] . 
x,(0 

is also written as 
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(14) 

where 

(15) 9 * = 2jni (fj Cos9 ji -ftji) 

Further, all multipaths from the mobile arrive at the base station array 
5 uniformly within ± A around the mean angle of the arrival Q. . 

Further, assuming the signals are narrowband, the complex baseband 

signal vector can be written as 

x,(0 « s/t - r 0 )(y L YidjyRje^ 0 ) 

(16) ^'=> ) 

Further letting 

10 ■yW = y 1 H0j)Rjie iv/t) 

(17) J ^^/=i J 

be an Array Response vector for L number of coherent multipaths 
corresponding to /th source. 

Use of the conventional algorithms like MUSIC leads to the estimation of 
the same direction 

15 0 j9 since all the L multipaths corresponding to that source are coherent. 

Hence 0 jt may be replaced by 0j . 



In such a case the above equation can be approximated as 



(18) 

For J-users the signal vector at the array is given by 

20 x(o = x x (0 + x 2 (o + + x y (o = V J Sj (ty&j + n(0 

(19) ^—^7=1 

n(t) is the noise vector (column), noise is assumed to be spatially and 
temporally white. 

(20) R = £[y(0y(0"] 
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R is a correlation matrix. It is obtained from the outer-products of the 
received signal+noise vectors x(t) and averaged over several snapshots. 
Averaging is denoted by E. 
H denotes the 'Hermitian(transpose and conjugate). 
A is the matrix of size NxJ (N-antennas, J-sources); 
a 2 is the noise variance, 
I is the NxN identity matrix and 

(21) 0> = E[s(t)s(t) H ] is the JxJ source covariance matrix. 

s(t) is baseband signal powers lying along the diagonal elements of O. 
It's size is JxJ, where J is the number of sources. Noise covariance is 



So, letting the eigenvalues of R, be arranged in descending order, and 
denoted by A, > X 2 >,... > X N , as well as letting u, ,...,u^ be the corresponding 

eigenvectors, where the eigenvalues X n are given by 



U s = [u, v ..,u y ] may be considered as the signal subspace, and 
U w =[u y+iv ..,u^] may be considered as the noise subspace. 
V n \(0j) = 0 for j = , and elsewhere U n v(0 y ) ± 0 . 

The embodiment assumes that the subspace spanned by the estimated 
eigenvectors U s =[u 1 ,...,u J ] corresponds to the space spanned by the J-true 

source directions. It should be noted that y'-th eigenvector corresponds to in- 
coherent multipaths due to the y'-th source. Estimation of the y-th source direction 
6 j can be done by any well known method like MUSIC. The embodiment 



given by 



(22) 



E[n(t)n H (t)] = a 2 n l 



(23) 




v n +o 2 forn = 1,...K 
a 2 for n = K +l,...N 
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searches within the direction range 0,. ±A , centered on the determined DOA of 
they-th source 0 y . , to determine the L multipaths for >th source. 

Thus, signal vector \j{t) can be considered as 

x,-(O = [v(0„) y(0 J2 ) ... v(^)]s}, 

(24) where s jt = R^e' 2 "^ cos V/'*> s . (, - tfl ) 

(25) x . (0 € span{y{6 n ),...v(0 yi )} . 

Thus, as illustrated in Fig. 3, the process for searching for the set 

0 ; = (Pji > d j2> •■•> d ji) where d ji e (Oj ± A )J = 1,2,.. .,L , starts with the process first 
selecting a trial direction set 0 = (d Jki ,d jki ,...,d Jki )for L multipaths of the yth 

source, block 302. 

The processing of finding 0 = 0 ; , that is, 

(@jk, >6jk 2 >-> e jk L ) = (0/1.0/2 >->9jl ) is explained as follows. 
First, the process selects randomly a set (0 jki ,d jki ,...,6 jk ) 

Second, the process forms the following matrix 

T>(0) = [v(e jki ) v(0 jk2 ) ... y{9 jkL ) u,] where 9 jki s(dj±A) 

Third, the process computes the following relationships 

(26) r u =\\v(d jki )\ 2 , and 

(27) q, =_Li^, block 304. 

||. || 2 is the 2-norm of the vector. 

Similarly, for following coefficients and vectors, the process computes 

(28) r a =qfy(e jki ); \<i<l-l; l = 2,...,L + \ 

(29) r„ =|| v(^, ( )-^/- f/ q,. || 2 , / = 2,...,Z + l,and 
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/-i 

(30) q, =[v(^ j )-SA} / q,]/r // , / = 2,...,L + 1 , block 306. The computation 
is performed until l=L+1, that is obtaining r ul Ui . 

Then, at block 308, the process computes the function 

(31)5(0,) = -!— . 

5 If the result of B(@j) yields a "new" set of peak values, the trial direction 

set 0 is set as the directions of arrival of the L strong multipaths. 

At block 310, a determination is made whether the process is to be 
repeated. The number of trials to repeat may be predetermined, for operating 
efficiency, or based on a predetermined threshold of diminishing marginal 
10 improvements. 

If another trial set is to be evaluated, the process returns to block 302. 

Eventually, the criteria to terminate estimation process is met, and the 
process terminates. The then current estimates of the L multipath directions 

A /V /\ 

&j = {0 jXy ...,6 jL } are used to obtain the combined signal z(0(See Figure.1) as 

L N 

15 (32) =ZZW where ™nl= e 

1=1 m\ 

In other words, after estimation of the L-multipath directional set in DOA 
Estimation Unit 112 fory'-th source, the received signals at the antenna elements 
are appropriately weighted as given in the equation (32). 

This process is repeated for all J signal sources received at the antenna 
20 elements with each source signal weighted by coefficients corresponding to L- 
multipaths. 

Referring now to Fig. 4, a block diagram illustrating an example computer 
system suitable for hosting a software implementation of the DOA Estimation unit 
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112 and Beamforming unit 114, is shown. As illustrated, computing device 400 
includes one or more processors 402, system memory 404, mass storage 
devices 406, other I/O devices 408 and network communication interface 410, 
coupled to each other via system bus 412 as shown. 
5 Processor 402 is employed to execute software implementations of DOA 

Estimation 112 and/or Beamforming 114. Processor 402 may be any one of a 
number of processors known in the art or to be designed. Examples of suitable 
processors include but are not limited microprocessors available from Intel Corp 
of Santa Clara, CA. 

1 0 System memory 404 is employed to store working copies of Estimation 

112 and/or Beamforming 114 and operating system services. System memory 
404 may be Dynamic Random Access Memory (DRAM), Synchronous DRAM 
(SDRAM) or other memory devices of the like. 

Mass storage devices 406 are employed to persistently store data, 

15 including e.g. a persistent copy of Estimation 112 and/or Beamforming 114. 

Examples of mass storage devices 406 include but are not limited to hard disks, 
CDROM, DVDROM, and so forth. 

Other I/O devices 408 are employed to facilitate other aspects of 
input/output. Examples of other I/O devices 408 include but are not limited to 

20 keypads, cursor control, video display and so forth. 

Network communication interface 410 is employed to facilitate network 
communication with other devices. Network communication interface 410 may 
be wired based or wireless. In various embodiments, network communication 
interface 410 may support any one of a wide range of networking protocols. 
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In alternate embodiments, computing system may be a multi-processor 
systems, a cluster of networked computers, including an array of massively 
parallel computing nodes. 

Accordingly, various novel methods and apparatus for determining the 
DOA of J signals, and/or determining the DOA of L multipaths of J signals, the J 
signals being wireless transmitted by their sources. 

While the present invention has been described in terms of the foregoing 
embodiments, those skilled in the art will recognize that the invention is not 
limited to the embodiments described. Other embodiments may be practiced 
with modification and alteration within the spirit and scope of the appended 
claims. Accordingly, the description is to be regarded as illustrative instead of 
restrictive. 
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